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We establish an information-theoretic approach for quantitatively characterizing the Non- 
Markovianity of open quantum processes. Here, the quantum Fisher information (QFI) flow provides 
a measure to statistically distinguish Markovian and non-Markovian processes. A basic relation be- 
tween the QFI flow and non-Markovianity is unveiled for quantum dynamics of open systems. For 
a class of time-local master equations, the exactly-analytic solution shows that for each fixed time 
the QFI flow is decomposed into additive sub-flows according to different dissipative channels. 
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I. INTRODUCTION 

Any system in the realistic world is open since it in- 
evitably interacts with its environment. The time evolu- 
tions of open systems are simply classified into Marko- 
vian and non-Markovian ones according to the ways to 
lose energy or information In most situations, Marko- 
vian process uniquely determines its final steady state as 
an thermal equilibrium [2j, which is independent of its 
initial one. In this sense a Markovian process is essen- 
tially an information erasure process, thus tends to con- 
tinuously reduce the distinguishability between any two 
initial states Q. 

However, Markovian description for an open quantum 
system is only an approximation to most of realistic pro- 
cesses, which are of non-Markovian. With many recent 
investigations about non-Markovian dynamics by mak- 
ing use of various analytical approaches and numerical 
simulations, a computable measure of "Markovianity" 
for quantum channels was introduced in Ref. [4]. Most 
recently, it was also recognized that difference between 
them can be measured through the continuous increment 
of the state distinguishability Q . Then this increment is 
intuitively interpreted as the revival of information flow 
between the bath and the system though there no quan- 
titative information measure is utilized. Based on this 
measure of the non-Markovianity, a method for direct 
measurement of the non-Markovian character was pro- 
posed [26|. Another approach based on entanglement is 
proposed in Ref. [5|. In this paper, the quantum Fisher 
information (QFI) flow are introduced to directly char- 
acterize the non-Markovianity of the quantum dynamics 
of open systems. 

Actually, in the system-plus-bath approach for open 
quantum systems, the effective dynamics of the reduced 
density matrix p is induced by tracing over environ- 
ment [1]. The simplest reduced dynamics is the quan- 
tum Markovian process described by dynamical semi- 



groups [6j. There, the reduced density matrix p at time 
t + dt is determined completely by the one at time t. 
Contrarily, the general reduced dynamics may be of non- 
Markovian when the surrounding environment may re- 
tain a memory of the information about states at earlier 
times, and transfer it back to the system to affect its evo- 
lution. In this sense the Markovian process only happens 
when the environmental correlation time is relatively 
short so that memory effects can be neglected. These 
memory-based considerations for the Markovian approx- 
imation also mean that the information-theoretical char- 
acterizing of the non-Markovianity is a quite natural fash- 
ion. However, it is still an open question that how to treat 
the information flow in open quantum systems based on 
a solid information-theoretic foundation. 

In this paper, we establish such a framework by adopt- 
ing the QFI flow as the quantitative measure for the in- 
formation flow. The QFI characterizes the statistical dis- 
tinguishability of reduced density matrix [7|, [8| . An in- 
tuitive picture of the memory effect of a non-Markovian 
behavior then immediately follows from the dynamic re- 
turn of the QFI, which is depicted by the inward QFI 
flow. For a class of the non-Markovian master equations 
in time- local forms, we exactly calculate the information 
flows. The analytic results show that the total QFI flow 
can be decomposed into the split contributions from dif- 
ferent dissipative channels for each fixed time. On the 
other hand, the QFI plays an essential role in quantum 
metrology [9[, where the highest precision of estimating 
an unknown parameter we may achieve is related to in- 
verse of the QFI. We point out this QFI flow approach is 
feasible to work for understanding the problems of quan- 
tum metrology. 



II. QUANTUM FISHER INFORMATION IN 
NON-MARKOVIAN DYNAMICS 



We consider the quantum processes described by the 
following time-local master equation 0, [Toj 
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d_ 

dt 



p(t) = K(t)p(t), 



(1) 
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where JC(t) is a super-operator acting on the reduced den- 
sity matrix p(t) as [1QI-KL2] 



)C(t)p=-i[H,p]+J2' 



A iP A] 



(2) 



with H (t) the Hermitian Hamiltonian for the open quan- 
tum system without the couplings to the bath. {•,•} 
denotes the anti-commutator. If all 7^ and Ai are 
time independent, and all ji are positive, equation (j2j) 
is the conventional master equation of the Lindblad 
form [11], which describes the conventional Markovian 
process. However, by making use of a variety of meth- 
ods, such as the time-convolutionless projection operator 
technique [l3) , Feynman- Vervon influence functional the- 
ory [14j and some others [15j, the parameters 7^ = ji(t) 
and Ai = Ai(t) in the time- local master equation may 
explicitly depend on time, and 7^ even may become nega- 
tive sometimes. Obviously, the non-Markovian character 
resides in these time-dependent coefficients. 

Taking some real number in the reduced density ma- 
trix p(Q- 1) as the inference parameter, we write down the 
QFI [H 



T(6;t) :=TY[L 2 (6;t)p(6;t)} 



(3) 



where L(0;t) is the so-called symmetric logarithmic 
derivative (SLD), which are Hermitian operators deter- 
mined by [16j 



^L p {9;t)= l -[L{9;t)p{e ] t) 



■p(6;t)L(6;t)}. (4) 



An important essential feature of the QFI is that we can 
obtain the achievable lower bound of the mean-square 
error of unbiased estimators for the parameter through 
the quantum Cramer-Rao (QCR) theorem 



Var(0;£) > 



1 



MF(6; t) ' 



(5) 



where M represents the times of the independent mea- 
surements [16] . From the QCR theorem, we can see that 
the QFI is indeed a measure of a certain kind of infor- 
mation with respect to the precision of estimating the 
inference parameter. The relations between the QFI and 
the statistical distinguishability of p(0; t) and its nei ghb or 
has been pointed out in some previous works [7|, [8|, [l6| . 

Flow of the QFI and its decomposition. — Here we use 
the QFI to characterize the non-Markovianity of the open 
quantum system by introducing the QFI flow, which is 
defined as the change rate X := dT jdt of the QFI. As 
a central result in this paper, a proposition about the 
decomposition of the QFI flow is given as follows: 

Proposition. For an open quantum system described by 
the time local master equation the QFI flow X = 
^2iXi is explicitly written as a sum of subflows Xi = jiJi 
with 



J := 



-Tr{p[L,Ai]i[L,Ai]}<0. 



(6) 



Proof From the differential of Eq. (j4]) with respect to 
time t, we have 



dtdep{9) = \ 



Lp + Lp + pL + pL 



It gives 



Tr 



pLL + pLL 



Tr[2Ld t d e p(6)] - Tr [2pL 2 



(7) 



(8) 



From the differential of Eq. (|3]) with respect to time t, 
we obtain the QFI flow as 



Tr 



(9) 



where the operator C := L(2d/d0 — L) is defined. By 
using the concrete expression of the master equation (j2j), 
we split the QFI flow into those individuals corresponding 
to the different dissipative channels as X = Tr [C}C(t)p(t)] 

or X = E; 7* Tr [£(^Mb - \£>i A \ A iiP}\- After some 
algebra, we get the decomposition X = ^iJi, where Ji 
is just given in Eq. (j6j). It finally proves the proposition. 

The above proposition and its proof contain rich impli- 
cations in physics. Firstly, the decomposition of the QFI 
flow corresponding to the different dissipative channels is 
due to the linearity of QFI flow equation (J9j) with respect 
to dp jdt and the concrete form of the time- local master 
equation (j2j). This is not a simple decomposition since 
each subflow depends on the whole SLD L{6]t), mean- 
while, L(0; t) is deduced from p(9; t), whose evolution de- 
pends on every dissipative channel and the unitary part 
of the master equation. So this kind of decomposition 
does not mean different dissipative channels are separa- 
ble to influence the change of the QFI for a period of 
time. However, for each fixed time t > 0, the QFI flow at 
the present moment are decomposed into the split contri- 
butions from different dissipative channels. In this sense, 
we interpret Xi(t) = 7i(£)j7i(£) as a subflow of the QFI 
at time t caused by the dissipative channel described by 
Ai(t) and ji(t). The magnitude of the QFI subflow is 
determined by a state-independent factor 7^ and a state- 
dependent factor Ji. 

Secondly, one of the advantages of such decomposition 
comes from the link between the direction of each QFI 
subflow Xi and the sign of the decay rate 7^. Because 
Ji is non-positive, we conclude that a negative ji(t) im- 
plies an inward QFI subflow (Xi > 0), except the triv- 
ial case of Ji = 0. The temporary appearance of neg- 
ative decay rates is already considered as the essential 
feature of the non-Markovian behaviors [l7j], here this 
is justified through the return of the QFI. For the case 
that all ji(t) are positive, the master equation ([2]) de- 
scribes a so-called time-dependent Markovian quantum 
process BEES El, in which cases, X always decreases. 
If the total QFI flow X(t) is positive at time £, it signifies 
at least one of ji(t) is negative. In such cases, the QFI 
flows back to the open system and the non-Markovian 
behavior emerges. 
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Actually, like the trace distance used in Ref. the 
dynamical return of the QFI is linked to the divisibility 
property of the dynamical map of quantum processes. If 
the master equation is of the form (|2j), the corresponding 
dynamical map is infinitely divisible provided that all 
ji are positive [20j. In such cases, for arbitrary time 
t > 0, the dynamical map from time t to t + dt is a 
completely positive and trace-preserving map. Thus the 
QFI decreases during this time interval since the QFI 
is monotonic with respect to a completely positive and 
trace-preserving map [21] . 

Thirdly, there should be some restriction on the evo- 
lution of the QFI. It is seen from the above proof of 
the proposition that the coherent part of Eq. (j2j), i.e. 
—i[H(t),p(t)], does not contribute to the total QFI 
flow directly. This observation directly leads to the no- 
cloning theorem in quantum information, which states 
that we can not use unitary operations to evolve the 
states \tb(Q)) |0) into \^{0)) \^{0)) as a quantum 
copy [22]. This is because the QFI of the target states is 
twice as the one of the source states, due to the additivity 
of the QFI for the product states. Besides, if the total 
system (system plus environment) is assumed closed and 
the QFI is only distributed in the system initially, then 
the QFI of the reduced density matrix during evolution 
should be not greater than the one at the initial time, for 
the invariance of the QFI under unitary evolution and 
the non-increasing of the QFI under partial trace opera- 
tion. This restriction should be reflected in the QFI flow 
obtained from a proper master equation. 
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FIG. 1: Estimation of parameter <j> in an unitary operation. 
After the phase gate operation, the system interacts with a 
reservoir. The precision of the estimation is impacted by char- 
acteristics of both the reservoir and the interaction. 



III. TWO-LEVEL SYSTEM 

Now we use an example of two-level system (qubit) to 
explicitly illustrate our discovery about the intrinsic re- 
lation between the QFI flow and the non-Markovianity 
of the open quantum system and its impaction to pa- 
rameter estimation. In the quantum metrology context, 
the QFI gives a theoretical-achievable limit on the preci- 
sion when estimating an unknown parameter, according 
to the QCR theorem ([5]). To estimate the parameter as 
precisely as possible, we should optimize input states to 
maximize the QFI, and then optimize measurements to 



achieve the Cramer- Rao bound [9]. However, due to the 
interaction with environment, the QFI will change and 
affect the precision of the parameter estimation. 

Here, the QFI-based parameter is assumed to be in- 
duced by a single-qubit phase gate := \g)(g\ + 
exp(z0)|e)(e| acting on the qubit, where <\> = is some 
inference parameters (see Fig. 1). To estimate the un- 
known parameter <p as precisely as possible, the optimal 
input state may be chosen as l^opt) = (\g) + l e ))/V% 
which maximizes the QFI of the output state E/^l^opt)? 
see Ref. [9[. In the following model, after the phase gate 
operation and before the measurement performed, the 
qubit is assumed as an atom coupled to a reservoir con- 
sisting of harmonic oscillators in the vacuum. The total 
Hamiltonian of this typical model p], [23|, HH reads 



H = u cr+cr- + ^ uj k b\b k + {a + B + (J-B^) (10) 
k 

with B = J2k9kbk, where uo denotes the transition fre- 
quency of the atom with ground and excited states \g) 
and |e), and a± the raising and lowering operators of 
atom; b\ and b k are respectively the creating and anni- 
hilation operators of the bath mode of frequencies uo k . 
g k denote the coupling constants. Then we consider 
Lorentzian spectral density J(u) = XW 2 /{tt[(ujo — uo) 2 + 
A 2 ]}, where W is the transition strength, and A defines 
the spectral width of the coupling, which is related to the 
reservoir correlation time scale tb by tb = A -1 p], |23| . 
The Lorentzian spectral density describe the reservoir 
composed of lossy cavity, see Ref. [1]. The time- local 
master equation of the form ([2]) can be obtained exactly 
as follows [l| 

^Ps(*) = lit) (*-p s (t)a + - l{ ff+ (T_,ps(t)}) (11) 

where j(t) = —2h(t)/h(t) with a crucial characteristic 
function 



h(t) = 



e- A£ / 2 [cosh(f) 
e -At/2 [cos(f)- 



f ^sinh(f)], 
Asin(f)], 



(12) 

where d = ^\X 2 -4W 2 |. 

Taking the initial state U(p\ip op t), the reduced den- 
sity matrix of the atom obeys the master equa- 
tion (fTTj) . Its solution is p s (t) = (I + B • a) /2, 
where B = (h(t) cos 0, — h(t) sin 0, h(t) 2 — 1) and a = 
(ctx^Oy^Gz). In order to calculate the QFI flow, we 
first diagonalize this reduced density matrix as ps(t) = 
^ZiPi(t)\^i(t))(^i(t)\. In this diagonal representation, 
the SLD with matrix elements = 2{ r ijji\d ( f ) ps\ r ^j) / {pi~\- 
Pj) is obtained explicitly as 

L(t) = ih(t) [\Mt))(Mt)\ - \Mt))(Mt)\] ■ (13) 

Further, we have J = -Tr(p [L, cr.] 1 [L,<t_]) = -h(t) 2 . 
Then the exact solution for the QFI flow 



l${t) = l(t)J(t) = 2h(t)h(t). 



(14) 
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is obtained, which leads to T§ = h(t) 2 . 

(a) (b) 
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FIG. 2: Two- level atom coupled to reservoir with Lorentzian 
spectral density: (a) QFI flow as a function of rescaled time, 
plotted in the weak coupling regime (W = 0.3A); (b) 7 as a 
function of rescaled time, W = 0.3A; (c) QFI flow as a func- 
tion of rescaled time, plotted in the strong coupling regime 
(W — 3A); (d) 7 as a function of rescaled time, W — 3A. 

Therefore, the characteristic of the QFI flow is deter- 
mined by the function h(t), which has two very differ- 
ent kinds of behaviors. The corresponding properties of 
the QFI flow are shown in Fig. ((2]). In the weak cou- 
pling regime (W < A/2), the function j(t) is always pos- 
itive, thus the QFI is always lost during the time evolu- 
tion of the open system. In the strong coupling regime 
(W > A/2), the function j(t) takes on negative values 
within certain intervals of time, see Fig. [2] (d), which 
displays the non-Markovianity. Obviously in these time 
intervals, the QFI flow is inward. It is remarkable that 
although j(t) diverges at certain time, the QFI flow does 
not, see Figs. [2] (c) and (d). This is because the QFI 



flow is determined by two factors, j(t) and J if). 

IV. CONCLUSION 

In summary, based on the QFI flow, we have pro- 
posed an information-theoretical approach for character- 
izing the time-dependent memory effect of environment 
for its surrounding open quantum systems. In this ap- 
proach the Markovian process is considered as a QFI era- 
sure process, and the return of the QFI, i.e. an inward 
QFI flow, clearly signatures the non- Markovian process. 
Using the time- local master equations, we have showed 
that for each fixed time t > the QFI flow is decom- 
posable according to different dissipative channels, and 
the direction of each sub-flow is determined by the sign 
of decay rates. With this decomposition form, the rela- 
tionship between the temporary appearance of negative 
decay rate and the non-Markovian characteristic is jus- 
tified. Although in the present work, the analysis of the 
QFI flow is based on a time-local master equation, the 
concept of the QFI flow may be still available in more 
general cases. 

The present approach is associated with the current 
development of quantum metrology, which concerns on 
finding an optimal fashion to make high-resolution and 
highly sensitive measurement of physical parameters 0. 
Due to the interaction with environment in experiments, 
like the photon losses in the optical interferometry or the 
presence of quantum noise [25], the QFI will change and 
affect the precision of the parameter estimation. There- 
fore, it is worthy to study the dynamical evolution of the 
QFI in the context of quantum metrology, especially for 
non-Markovian processes. 
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